Nonlinear Optical Spectroscopy of Photonic Metamaterials 
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We have obtained spectra of second-harmonic generation, third harmonic generation, and four- 
wave mixing from a fishnet metamaterial around its magnetic resonance. The resonant behaviors are 
distinctly different from those for ordinary materials. They result from the fact that the resonance 
is plasmonic, and its enhancement appears through the local field in the nanostructure. 

PACS numbers: 



Optical metamaterials with nanoscale metal building 
blocks have been studied extensively in recent years 
Q. With each metal unit much smaller than optical 
wavelength, they can be viewed as continuous media and 
the metal units act as 'artificial molecules'. The optical 
properties of metamaterials can be engineered by proper 
design of the 'artificial molecules' and they can exhibit 
unusual behavior nonexistent in nature, such as nega- 
tive refractive indices. While linear optical properties of 
metamaterials have been well investigated [2j]- [H, non- 
linear optical properties began to attract interest only re- 
cently 0- Such interest stems from possible strong 
enhancement of nonlinear response from plasmon reso- 
nances of the metal nanostructures together with spec- 
tral tunability offered by design of 'artificial molecules'. 
A crucial aspect in understanding nonlinear optical prop- 
erties of metamaterials is their spectral responses, which 
have not yet been reported. 

In this letter, we present the first spectroscopic study 
of second harmonic generation (SHG), third harmonic 
generation (THG) and four-wave mixing from a meta- 
material comprising a monolayer of "fishnet" struc- 
ture [2fj| . [2lj . It was designed to have negative refrac- 
tive index in the near-IR range, with a magnetic reso- 
nance around 1.55 /im 22]. The spectra of SHG and 
THG with the fundamental input scanned over the mag- 
netic resonance were obtained with different fundamental 
and harmonic polarizations. Resonant enhancement were 
clearly observed. Interestingly, the observed resonances 
are much sharper than that in linear absorption. This 
is distinctly different from typical molecular cases, where 
resonant excitation at the fundamental wavelength yields 
the same resonance spectrum in linear and nonlinear re- 
sponses. Such difference originates from the fact that, 
unlike molecular resonances, the plasmon resonances in 
metal nanostructures are collective oscillations and their 
resonance enhancement appears through the local field 
effect in the nonlinear processes. 

The measurements were carried out on a "fishnet" 
metamaterial composed of two silver sheets with hole 
arrays separated by a Si02 layer. It was fabricated us- 



ing combination of nanoimprint lithography (NIL) and 
electron-beam lithography (EBL) [23[. The SEM image 
of the structure is shown in Fig. la and the structural 
configuration of the broad wire and its dimensions in Fig. 
lb. Linear optical response of this metamaterial exhibits 
a magnetic resonance at 1.55 /im when the magnetic field 
component of the input wave threads the loop formed by 
linking the broad metal wires of the two layers, as indi- 
cated in Fig. lb by the black arrows [3], [25[. Corre- 
spondingly, the effective refractive index is negative in the 
wavelength range of 1.45 to 1.6 /jm. The measured lin- 
ear transmittance and reflectance spectra are presented 
in Fig. lc. 
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FIG. 1: a) SEM image of the fishnet structure; b) Schematic 
structure of the broad wire of "fishnet" Ag/SiC>2/Ag struc- 
ture, c) Linear transmittance (open dots) and reflectance 
(solid dots) spectra. 

The nonlinear optical spectroscopy of the fishnet struc- 
ture was performed using the tunable output from an op- 
tical parametric system, pumped by the third-harmonic 
of a picosecond YAG:Nd laser. The tunable IR beam was 
incident at 30° on the sample, and the reflected SHG, 
THG and four wave mixing signals were detected by a 
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photo-multiplier and gated electronics system after spec- 
tral filtering. 
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FIG. 2: SHG spectra of the fishnet structure for different 
polarization combinations: (a) p-in p-out, (b) p-in s-out, (c) 
s-in s-out, (d) s-in p-out. Solid dots and open dots are for 
beam geometry with the incident magnetic field component 
along and perpendicular to the broad Ag wire, respectively. 

The SHG spectra, obtained with the incident plane 
along (solid dots) and perpendicular (open dots) to the 
broad Ag wires of the fishnet are shown in Fig. 2. Polar- 
ization combinations employed are indicated in the fig- 
ure (PS, for example, denotes P and S polarizations for 
the fundamental input and SH output, respectively). To 
correct for the wavelength dependent incident laser inten- 
sity, the SHG signals at each wavelength were normalized 
to that of a smooth silver film with the PP polarization 
combination. Two features are clear in the SHG spec- 
tra. First, the SHG spectra display a resonance at 1.55 
fim whenever the input polarization has a magnetic-field 
component along the broad wires of the fishnet (i.e., P- or 
S-polarized when the incident plane is perpendicular to 
or along the broad wires, respectively) and the magnetic 
resonance is excited. Otherwise, the SHG spectra are 
featureless with only non-resonant contribution. Second, 
the PS and SS polarizations yield weaker SHG signals 
with a much lower nonresonant background compared to 
PP and SP. Presumably this is the result of symmetry. 
As in the case of thin Ag films, SHG from a perfect fish- 
net structure, with the incident plane coinciding with a 
mirror plane, is strictly forbidden for PS and SS polar- 
izations, but allowed for PP and SP. That SHG with SS 
and PS is actually observable is an indication that the 
fishnet sample is not ideally symmetric. 

In Fig. 3, the THG spectra are shown for different 
sample orientations and polarization combinations using 
the same notations as those for SHG in Fig. 2. Strong en- 



hancement of THG is observed when and only when the 
magnetic resonance is excited by the fundamental beam. 
The stronger signal for PP and SS compared with PS and 
SP polarization can also be understood from symmetry 
argument: THG with PS and SP is forbidden in a perfect 
fishnet structure when the incident plane coincides with 
a mirror plane. Compared with the SHG spectra, the 
non-resonant contribution is smaller in the THG signal 
and the resonance more pronounced. 

We compare in Fig. 4 the magnetic resonant features in 
the PP spectra of linear absorption, SHG and THG. The 
resonance of THG is clearly narrower than that of SHG, 
which is in turn sharper than that of linear absorption. 
This is in striking contrast with molecular resonances, 
where resonance lineshapes of harmonic generation spec- 
tra are similar to that of linear absorption. It indicates 
that resonances of "artificial molecules" in mctamatcri- 
als are characteristically different from those of natural 
molecules. The difference arises because plasmon reso- 
nances of metal nanostructures are intrinsically collec- 
tive in nature instead of local as in molecular transitions. 
Such distinctions do not show up in linear optical prop- 
erties, but become obvious in nonlinear optical spectra 
where a resonant input field participates multiple times 
in nonlinear processes. To further establish the above- 
mentioned characteristic resonant behavior in nonlinear 
responses of metamaterials, we measured outputs of four 
wave mixing at 3u>i, 2loi+lu2 and lui+2(V2 with u>2 fixed 
(at wavelength of 1064 nm) and u% scanned across the 
magnetic resonance. The observed spectra for PP polar- 
ization are displayed in Fig. 4b with signals at 2lui+lu2 



» PP 
:■ ■Hi-* 1 * 


» k 


A 


d) SP 



1.56 1.60 1.64 152 1.56 1.60 
Fundamental wavelength (Tin) 



FIG. 3: THG spectra of the fishnet structure for different 
polarization combinations: (a) p-in p-out, (b) p-in s-out, (c) 
s-in s-out, (d) s-in p-out. Solid dots and open dots are for 
beam geometry with the incident magnetic field component 
along and perpendicular to the broad Ag wire, respectively. 
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and toi+2u2 scaled up by 6.6 and 43, respectively. It is 
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FIG. 4: a) Comparison of SHG (triangle) and THG (open 
dots) spectra in pp polarization combinations with the linear 
absorption spectrum (solid dots), b) Spectra of four- wave 
processes at 3o;i (open dots), 2ui+lj2 (triangles) and wi+2w2 
(solid dots). The signal strengths at 2^1+^2 and uji+2ui2 are 
scaled up by 6.6 and 43, respectively, for ease of comparison. 

obvious that each additional cjj component in the input 
leads to extra enhancement of the output at resonance, 
and the corresponding resonance peak becomes increas- 
ingly sharper. 

For better understanding of the observed results, we 
realize that the real source of nth harmonic generation 
in metamaterials is the induced effective electric dipole 
on each nanostructural unit: 

P [n) = [ L(f,nw) : x (n) (r) : [E loc (r,u)] n dV, (1) 
Jv 

where the integration is over the volume of the unit, 
X {n} is the local rath-order nonlinear susceptibility and 
Ei oc (r, u) = L(r, u>) : Eq(uj) is the local field with L(r, u>) 
being the local field correction factor and Eq(uS) the in- 
put field. For simplicity, we have neglected the contribu- 
tion of magnetic-dipole and electric-quadrupole to P(") 



in Eq.(l). The harmonic output with polarization along 
n is given by oc \h ■ p(")| 2 . The local field here 

can be decomposed into two components, one associated 
with the magnetic resonance and the other not. The res- 
onant component is relatively enhanced when approaches 
resonance. Thus we can write the local field correction 
factor as L(r,uj) = A(r, u>) + B(r,u>)/D(uj) with D ap- 
proximated by D = u> — luq + iT describing the magnetic 
plasmon resonance. 

The symmetry requirement for harmonic generation 
is naturally incorporated in the volume integration of 
Eq.(l), which, for example, vanishes for symmetry for- 
bidden processes in a perfect fishnet structure. The in- 
tegral depends sensitively on the field distribution, and 
nonresonant and resonant terms respond differently to 
nanostructure change because of their different local field 
distributions: For a symmetry- forbidden harmonic gen- 
eration process, if symmetry-breaking modification of the 
nanostructure is not severe, the nonresonant signal is still 
expected to be small. At resonance, however, B(f,u>) be- 
ing different from A(f, ui) can have a spatial distribution 
emphasizing contribution from the symmetry-breaking 
part of the structure [26[, thus generating a relatively 
strong resonant harmonic output. This explains our ob- 
servation of resonant spectra with very weak nonresonant 
background for symmetry-forbidden harmonic generation 
processes in the fishnet structure presented in Figs. 2 and 
3. Numerical calculation on a realistic fishnet structure 
hopefully will quantify the different effects of symmetry 
breaking on resonant and nonresonant harmonic genera- 
tion. 

If the nonresonant part of L(r, ui) could be neglected, 
we would have, for the nth harmonic generation, p(") oc 
\D\- n and oc |L>|- 2 ". With the presence of the 

nonresonant A(f,oS) term in L(f,u>), the signal S^™' now 
has terms of |Z?|~ 2 ™ with m = 0, 1, n, and its resonant 
lineshape often appears broader than |Z?| -2 ™. However, 
the |D| -2 " term is always significant, making the THG 
(n = 3) spectrum sharper than that of SHG (n = 2); 
similar behavior is seen in comparing spectra of wave 
mixing at 3wi, 2u>i+uj2 and u)i+2u>2- 

In summary, we have measured spectra of SHG, 
THG and four-wave mixing from a fishnet metamate- 
rial around its magnetic resonance with different in- 
put/output polarization combinations. The results show 
that the resonant enhancement is much stronger in non- 
linear optical response than that in the linear case, and 
the more times the resonant input field participates in 
the mixing process, the sharper the resonant spectrum 
appears to be. This is because the resonance is plasmonic 
in nature and shows up in the local-field correction fac- 
tors in the nonlinear responses. Such resonant behavior 
is expected to appear in all nonlinear optical processes 
involving plasmon resonances in metamaterials. 
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